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The electric dipole moment of the neutron (nEDM) is a very sen- 
sitive probe for CP violation beyond the Standard Model of par- 
ticle physics PPJ2] • The most precise experiments today use Ram- 
sey's resonance method and ultra-cold neutrons (UCNs) 0@]. 
Further progress is presently limited by systematics [5j and the 
low density of UCNs available. Here we discuss an alternative ap- 
proach based on spin rotation in non-centrosymmetric crystals. 

The statistical sensitivity of any experiment to measure the nEDM 
is determined by the product Et\/~N \ where r is the duration of 
the neutron interaction with the electric field E and N the num- 
ber of the counted neutrons. New projects to measure the nEDM 
with UCNs aim to increase the UCN density and thus N by or- 
ders of magnitude (see [6j for a recent overview). In contrast, 
experiments with crystals exploit the electric field inside matter, 
which for some crystals can be by a few orders of magnitude 
higher than the electric field achievable in vacuum. 
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EDM experiments with absorbing crystals were pioneered by 
Shull and Nathans J7J. Their experiment was based on the in- 
terference of the electromagnetic amplitude with the imaginary 
part of the nuclear one. Abov with his colleagues [Sj were the 
first who paid attention to the presence of a spin dependent term 
due to the interference of nuclear and spin-orbit parts of the 
scattering amplitude in the interaction of neutrons with a non- 
centrosymmetric non-absorptive crystal. Spin-rotation in non- 
centrosymmetric crystals due to such interference effects as a way 
to search for a nEDM was first discussed by Forte [Sj. The cor- 
responding spin-rotation effect due to spin-orbit interaction was 
experimentally tested by Forte and Zeyen |10l| . The authors of 



PTTlfT2] have shown that the interference of the nuclear and the 
electromagnetic parts of the scattering amplitude leads to a con- 
stant strong electric field, acting on a neutron during all time 
of its movement in the non-centrosymmetric crystal. This field 
was measured in a Laue geometry diffraction experiment [12], in 
agreement with the calculated value. 

The spin rotation can be measured in Bragg [9 1fTUlfT3] and Laue 
[T^lfToTfTBlfTT] diffraction geometry. In this paper, we compare both 



geometries and present preliminary results of a test experiment in 
Bragg geometry. We show that the sensitivity of an optimized ex- 
periment in Bragg geometry can compete with the most sensitive 
published UCN nEDM measurements. 



2 Comparison of Laue and Bragg diffraction geometry 



A detailed recent study of a nEDM measurement in Laue geom- 
etry can be found in [18] . The main advantage of this scheme is 
the possibility to increase the time r of neutron passage through 
the crystal using Bragg angles 9b close to tt/2 [14]. In this way, 
times close to the time of neutron absorption r a ~ 1 ms were 
obtained in short crystals |T6l[T9] . However, in order to suppress 
systematic effects due to the Schwinger interaction (spin-orbit 
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coupling), the method relies on an effective depolarization of the 
neutron beam by Schwinger-rotating the two components of the 
neutron wave by ±7r/2. This fixes the thickness of the crystal to 



_ 7rm p c 



Eg is the electric field affecting the neutron for the exact Bragg 
condition for the crystallographic plane g (g is a reciprocal lattice 
vector) and m p the proton mass [T2] . Therefore the sensitivity of 



the method cannot be increased by using crystals with a higher 
electric field E g since the crystal length and thus the time r 
would need to be reduced. On the other hand, it is impossible 
to increase the sensitivity using Bragg angles extremely close to 
7r/2, because of a resulting decrease of the neutron count rate for 
such angles [IB] • In the experiment [IB] , also a strong sensitivity 
of the method to crystal deformations was found. 



The main advantage of the Bragg diffraction scheme [9,13j is that 
the electric field acting on the neutron depends on the deviation 
of the neutron trajectory from the Bragg condition. This allows 
us to control value and even sign of the electric field and makes 
new tests of systematic effects possible. On the other hand, the 
time r that the neutron spends in the crystal cannot be increased 
by using Bragg angles close to tt/2 as it depends on the total 
neutron velocity v and not on the velocity component parallel 
to the crystallographic planes as in the Laue case. However, this 
disadvantage can be obviated by increasing the crystal thickness, 
in principle. 



A first experiment in Bragg geometry was reported in [TU 



mea- 
suring the neutron spin rotation due to Schwinger interaction in 
quartz. The experimental value of the spin rotation angle was a 
few times less than the theoretical expectation [ID] . The main ex- 
perimental difficulty was to obtain monoenergetic neutrons with 
a well-defined small deviation from the Bragg condition. This was 
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solved by placing the monochromator in a strong magnetic field, 
but this field may cause systematic errors. 

Here we propose and use a very simple solution of the problem 
to obtain these monoenergetic neutrons. 



3 Two-crystal scheme 



Let's consider the symmetric Bragg diffraction case. A neutron 
falls on the crystal in a direction close to the Bragg one for the 
crystallographic plane g. The deviation from the exact Bragg 
condition is described by the parameter AE" g = Ek — Ek g , where 
Ek = h 2 k 2 /2m and Ek g = h 2 \k + g\ 2 /2m are the energies of a 
neutron in the states \k) and \k +g), respectively. 

In this case the neutron wave function inside the crystal in first 
order of perturbation theory can be written [2D] 

ij}(r) = e- ikr + a ■ e- i(fc+9)r , (2) 



where 



a 1^1 W , (3) 
Ek — Ek g AE g 



Here V g is the g-harmonic of the neutron-crystal interaction po- 
tential. For simplicity we consider the case a< 1, permitting to 
use perturbation theory. 

The electric field acting on the neutron in the crystal is [20] 
E = E g -a : (4) 



where E g is the interplanar electric field for the exact Bragg 
condition. A nonzero nEDM d n results in neutron spin rotation 
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by the angle 

2E-d n -L 

^EDM = Z , (5) 



where L is the length of the crystal and v± the component of 
the neutron velocity perpendicular to the crystallographic plane. 
Obviously, sign and value of the electric field (Hj) are determined 
by sign and value of the deviation from the exact Bragg condition 
AE g . This is illustrated in Fig. CD 

The presence of the electric field will lead to the appearance of a 
Schwinger magnetic field 

H s =[Ex v]/c : (6) 



where v is the neutron velocity. The corresponding spin rotation 
angle is 



^ = ~lwT' (7) 



Note that the Schwinger effect disappears for Bragg angles of tt/2 
as E || v in this case: 

2EuLv\\ 2EuL 
¥>s = ; " =^^cotfl B —►,(), 8 



where vn is the component of the neutron velocity parallel to the 
crystallographic plane. This can be used in the nEDM experiment 
to suppress systematic effects due to Schwinger interaction. 

In the experiment [21] the effect of neutron spin rotation in neu- 
tron optics for large (~ 10 3 — 10 4 Bragg widths) deviations from 
the exact Bragg condition have been observed. The measured 
effect has coincided with the theoretical expectation. 

In order to tune the electric field E, we have to select neutrons 
with a well-defined deviation AE g from the Bragg condition. For 
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this purpose, we use two separate crystals in parallel orientation 
(see Fig. [2J). By heating or cooling the second (small) crystal the 
interplanar distance Ad changes and, therefore, the energy of the 
reflected neutrons. The second crystal selects only the neutrons 
corresponding to its own Bragg condition from the whole beam 
passing through the first crystal. The deviation parameter and, 
accordingly, value and sign of the electric field having affected 
these neutrons in the first crystal depend directly on the temper- 
ature difference between the two crystals (see Fig. CD). 




-I ► 

TC/2 Bragg angle 



Fig. 1. Electric field for neutrons close to the Bragg condition (Bragg angle ir/2 
and A = 2do) in a non-centrosymmetric crystal: The electric field has a pole for 
neutrons exactly fulfilling the Bragg condition. By choosing the deviation from the 
Bragg condition, value and sign of the electric field can be selected. We use Bragg 
reflection in a second crystal with a different temperature To± AT to select neutrons 
with a certain deviation. 

The value of the Bragg width AAb for the (110) quartz plane 
(d = 2.45A) is AA B / A ~ 10~ 5 . To shift the reflected wavelength 
by one Bragg width, Ad/d should have the same value, corre- 
sponding to a temperature difference of AT ~ ±1 K (linear coef- 
ficient of thermal expansion for quartz £ = AL/L « 10~ 5 K _1 ). 
Note that a common variation of the two crystals' temperatures 
in itself does not influence the value of the electric field; only 
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the temperature difference between the two crystals needs to be 
controlled. 




Fig. 2. Two crystals are in parallel position. The crystal deviation parameter AE g 
for the neutrons reflected by the small crystal is determined by the temperature 
difference AT. 

4 Experimental setup 

A test experiment was carried out at the facility for particle 
physics with cold neutrons PF1B [22] of the Institut Laue-Langevin 
The scheme of the experiment is shown in Fig. [31 We used the 
(110) reflex of quartz. The length of the main crystal was 14 cm. 
As mentioned above the backscattering geometry (Bragg angle 
close to 7r/2) allows us to suppress effectively the Schwinger ef- 
fect. 

To measure the nEDM we should direct the initial neutron spin 
along the X axis (see Fig. |3J) and investigate the dependence of 
the Y component of the final polarization vector (i.e. the element 
Mxy of the polarization tensor) on the temperature difference 
between the two crystals. The main systematic error of this ex- 
periment is due to the residual Schwinger effect. It can be investi- 
gated by measuring the polarization component along the Z axis 
(elements Mxz and Myz of the polarization tensor). Therefore 
the measurement of a single projection of the final polarization 
vector only is not sufficient. We used the 3-D polarization anal- 
ysis device CRYOPAD [23] to measure the whole polarization 
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tensor. In first order the measured difference of the polarization 
tensor after reversing the electric field is 



AM = g n r 
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where To = (r+ + t_)/2, Ar = (r + — r_)/2, and r + and r_ are 
the times the neutrons stay in the crystal for the positive and the 
negative electric field, respectively. H % are the components of the 
residual magnetic field and H l s the components of the Schwinger 
magnetic field H§. g n = 2/i n /h = 1.8 • 10 4 G _1 s _1 is the neutron 
gyromagnetic ratio, #edm = Ed n / fi n the effective magnetic field 
corresponding to the electric field E and the nEDM d n . For E = 
1 • 10 8 V/cm and d n = lO" 26 e-cm, H EBM = 1.7 • 1CT 7 G. 

Monitor Beamstop 



Beamstop 




Fig. 3. Scheme of the experiment. The neutrons, coming from the right, pass a py- 
rolytic graphite monochromator (adjusted to 4.91 A), a super mirror polarizer, and 
a resonance spin flipper. The incoming nutator and a coil inside the Meissner cavity 
orient the incident neutron polarization vector in space before the main crystal. 
Correspondingly, a second coil and the outgoing nutator select the component of 
the outgoing neutron polarization that is analyzed by a 3 He spin filter. The ana- 
lyzer crystal back-reflects neutrons according to its Bragg condition. The pyrolytic 
graphite PG (002) with about 50 % reflectivity directs the reflected neutrons to the 
nPSD. It also unavoidably reflects a part of the beam after the spin filter which was 
used to monitor the beam intensity. 

Two main effects can simulate a nEDM: 
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(1) Schwinger effect. It results in nonzero components AM X z = 
-AM ZX « Hig n r and AM YZ = -AM ZY « #f# n r but 
does not contribute directly to the components AMxy = 
—AMyx related to the nEDM. However, false effects can 
arise from imperfections of the 3-D polarization analysis de- 
vice. These effects can be investigated by measuring the resid- 
ual Schwinger effect via AM X z or AM ZX - For example, if 
the angular accuracy of the 3-D analysis is 10~ 3 and the 
Schwinger effect is suppressed down to a level of 10~ 4 G 
(quite realistic numbers), then the residual false effect will 
be on the level 10~ 7 G, corresponding to d n ~ 6 • 10 e-cm. 

(2) Residual magnetic field. Effects due to nonzero magnetic field 
do not compensate in the difference tensor since the neutrons 
for both signs of the deviation parameter AE g have slightly 
different energies and thus spend slightly different times in 
the residual magnetic field: At/ 'tq = —A A/ A. In our case 
AA/A « 10~ 5 . Therefore the residual magnetic field in the 
polarization analysis device needs to be small: 

H*— « //EDM- (10) 
TO 

The main idea to control the contribution of the Schwinger effect 
to the nEDM matrix element AMxy is to measure the angular, 
i.e. the spatial distribution of the polarization tensor by a posi- 
tion sensitive neutron detector (nPSD) because the nEDM and 
Schwinger effects have different angular dependences for Bragg 
angles close to w/2. The nEDM effect is constant, see (ED, whereas 
the Schwinger one is proportional to cot(#B) ~ ( 7r /2 — 6*b), see (ED 
and would be visible in a spatial dependence of the components 
AM XY and AM YX - 

5 Preliminary results 

The measured dependence of the integral intensity of neutrons 
reflected by the second crystal on the temperature difference AT 
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Fig. 4. Dependence of registered intensity of neutrons on the temperature difference 
AT between the two quartz crystals. 

between the two quartz crystals is shown in Fig. HI The intensity 
minimum corresponds to equal interplanar distances in the two 
crystals. The observed difference of interplanar distances for the 
two crystals at AT = is probably due to impurities. 

In the previous measurements [21] we found the temperature dif- 
ferences corresponding to maximal electric fields acting on the 
neutron: AT = -2.0 K and AT+ = +0.4 K (arrows in Fig. 1). 
The corresponding electric field was determined from the slope 
of the spatial dependence of AM Z x (see [25] for the graph): 



This value coincides with the preliminary result obtained for the 
Bragg angle 6 B = 86° g3j. 

The analysis of the spatial dependence of matrix element AMyx 
showed no contribution from the Schwinger effect within the 
present statistics. Summarizing all data, we find the value for 
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the angle of neutron spin rotation related to the nEDM: 

A^edm = (0.9 ± 2.3) • 10- 4 . (12) 



From this value, we derive the preliminary result 

d n = (2.5 ± 6.5) • 10~ 24 e • cm. (13) 

The error is only statistical. However, from the preliminary anal- 
ysis of systematic errors we expect that the result is limited by 
statistics. The main contributions to the systematic uncertainty 
are due to the Schwinger effect, the residual magnetic field inside 
CRYOPAD, a small misalignment between the two crystals, and 
the curved exit window of CRYOPAD. 

6 Statistical sensitivity of a dedicated setup 

The experiment described above used existing equipment. The 
count rate in the "gray" position of the polarization measure- 
ments was only ~ 60 n/s and the statistical sensitivity to the 
nEDM was 1.6 • 10 -23 e-cm per day. The main limitations were: 

• Limited divergence acceptance. The beam divergence accepted 
by the installation was about ±2.5 • 10~ 3 only. However, the 
beam divergence after the polarizer is about ±1.5 • 10~ 2 . By 
installing neutron guides between the different elements, by 
shortening the installation, and by using a CRYOPAD with 
larger angular acceptance a factor of in statistical sensi- 
tivity can be gained. 

• Limited beam size. The used beam cross-section 5.5 cm 2 was 
limited by the size of the monochromator, the windows of CRY- 
OPAD, and the size of the crystals used. In a dedicated setup 
and using larger available crystals, a beam cross-section of 
6 x 12 cm 2 is feasible, yielding a factor of y/lS in statistical 
sensitivity. 
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• Limited crystal length. The interaction time of the neutrons 
with the electric field can be increased up to the absorption 
time constant in quartz, by using one large single crystal or 
a sequence of shorter crystals. We aim for a crystal length 
of 50 cm compared to 14 cm used in the experiment. Taking 
into account the absorption, this will yield a factor of 2.8 in 
statistical sensitivity. 

• Background, polarization analysis. Optimization of background 
and polarization analysis can yield a factor of 1.5 in statistical 
sensitivity. 

With these improvements, we expect an improvement of a factor 
of about 65 in statistical sensitivity, yielding a sensitivity of about 
2.5 • 10~ 25 e-cm per day. This is comparable to state-of-the-art 
nEDM measurements using UCNs. 

7 Conclusions 

We have measured the nEDM by crystal diffraction in Bragg 
geometry. Our preliminary result is d n = (2.5 ± 6.5) • 10~ 24 e • 
cm (only statistical error). The statistical sensitivity of the test 
experiment was 1.6 • 10~ 23 e-cm per day. It can be improved by a 
factor of 65 in a dedicated setup and can reach ~ 2.5 • 1CT 25 e-cm 
per day for available quartz crystals and neutron beam flux. The 
present scheme of the experiment allows to control systematic 
effects related to the Schwinger effect on-line with an accuracy 
better than the statistical one. The technical requirements (e.g. 
accuracy of 3-D polarization analysis, alignment of crystals) for a 
systematic uncertainty on the level of 10 -26 e-cm seem achievable 

121. 
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